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Title : Method and Apparatus for Determining Peripheral Breast 
Thickness 

FIELD OF THE INVENTION 

5 This invention relates in general to a method and apparatus for 

determining the thickness of a breast subjected to a mammogram, and more 
specifically relates to a method and apparatus for determining the thickness of 
a breast at its peripheral portion. 

BACKGROUND OF THE INVENTION 

10 In conventional mammography, a woman places her breast on a 

breast support plate of the mammography machine. A detector is typically 
mounted under the breast support plate. This detector is sensitive to x-rays. A 
breast compressor plate that is transparent to light and x-rays presses down 
against the top of the breast to flatten it and to prevent any movement of the 

15 breast during the mammogram. An x-ray source is then turned on to image 
the breast between the breast support plate and the breast compression 
plate. 

Mammograms provide clues that help to distinguish benign and 
malignant breast diseases. Radiologists look at both the static appearance of 

20 the breast, as well as changes in its structure, micro-classification, density 
and other characteristics. Breast density determined from the mammogram 
has been linked to increased link of breast cancer. Women with high 
mammographic densities (i.e., a high proportion of radiographically-opaque 
stroma and parenchyma) have been shown to be at an increased risk of 

25 breast cancer, when compared to a woman whose breasts are composed 
mainly of fatty or adipose tissue. Classification of radiological appearance of 
mammograms on the basis of the general distribution of parenchyma, stroma 
and fat, can yield very strong estimates of breast cancer risk. 

In the mammography field, various systems and methods have 
30 been developed to quantify breast density in terms of the fraction of the 
projected breast area that is occupied by radiographically dense tissue. These 



-2- 

methods suffer from at least two limitations. First, they do not use information 
about three-dimensional conformation of the breast. A simple area 
measurement may provide an erroneous measure of the actual amount of 
fibroglandular tissue in the breast. 
5 The computation of volumetric density in a compressed breast is 

based on both image data and knowledge of the thickness at each pixel. 
However, at the breast periphery, where the breast is not bounded by either 
the breast support plate or the breast compression plate, the thickness of the 
breast may not be known. However, this thickness is required to determine 
10 volumetric density of the compressed breast. 

Accordingly, a method and apparatus for determining the 
thickness of a breast at its periphery is desirable. 

SUMMARY OF THE INVENTION 

An object of one aspect of the present invention is to provide a 
15 method of generating a three-dimensional breast thickness object for a digital 

mammogram of a breast. 

In accordance with one aspect of the present invention, there is 
provided a method of generating a three-dimensional breast thickness object 
for a digital mammogram of a breast. The method comprises: 
20 (a) generating a phantom thickness object for transforming into 

the breast thickness object, the phantom thickness object being generated in 
a three-dimensional modeling means and being substantially breast-shaped; 

(b) determining a set of dimensions for the breast; and 

(c) transforming the phantom thickness object to conform to the 
25 set of dimensions to provide the three-dimensional breast thickness object in 

the three-dimensional modeling means. 

An object of a second aspect of the present invention is to 
provide a computer program product for use on a computer system for 
analyzing digital mammograms. 
30 In accordance with a second aspect of the present invention, 

there is provided a computer program product for use on a computer system 
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for analyzing digital mammograms. The computer program product 
comprises: 

(a) a recording medium; 

(b) phantom thickness object generation means recorded on the 
5 recording medium for instructing the computer system to generate the 

phantom thickness object; 

(c) data entry generation means recorded on the recording 
medium for instructing the computer system to upload a set of dimensions for 
the breast; and, 

-JO (d) transformation generation means recorded on the recording 

medium for instructing the computer system to transform the phantom 
thickness object to conform to the set of dimensions for the breast to provide 
the three-dimensional breast thickness object 

An object of a third aspect of the present invention is to provide 
15 a computer system for analyzing digital mammograms. 

In accordance with a third aspect of the present invention, there 
is provided a computer system for analyzing digital mammograms. The 
computer system comprises: 

(a) phantom thickness object generation means for generating 

20 the phantom thickness object; 

(b) data entry means for receiving a set of dimensions for a 

breast; and, 

(c) transformation means for transforming the phantom 
thickness object to conform to the set of dimensions for the breast to provide 

25 the three-dimensional breast thickness object. 
BRIEF DESCRIPTION OF THE DRAWINGS 

A detailed description of preferred aspects of the invention is 
provided herein below with reference to the following drawings, in which: 

Fig. 1, in a perspective view, illustrates a mammography 

30 machine; 
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Fig. 2, in a sectional view, illustrates a breast phantom image 
constructed of poly-methyl-methacrylate (PMMA); 

Fig. 3, in a perspective view, illustrates a three-dimensional 
triangle phantom image; 
5 Fig. 4, in a schematic view, illustrates a breast compressed 

during a mammogram; 

Fig. 5 is a graph of polynomial conversion functions obtained 
from the three-dimensional triangle phantom image of Fig. 3; 

Fig. 6 is a graph of a grey value histogram of a digital 
10 mammogram; 

Fig. 7 shows a phantom thickness map object; 

Fig. 8 shows a digital mammogram object; 

Fig. 9 illustrates the phantom thickness map object of Fig. 7 
including internal and external sets of landmarks; 
15 Fig. 10 illustrates the digital mammogram object of Fig. 8 

including internal and external sets of landmarks; 

Fig. 11 is a graph of a thickness profile for the phantom 
thickness map object of Fig. 7; 

Fig. 12 is a graph illustrating a thickness profile of the digital 
20 mammogram of Fig. 8; 

Fig. 13 is a flowchart illustrating a method of generating a 
phantom thickness map in accordance with an aspect of the invention; 

Fig. 14 is a flowchart illustrating a method of generating a breast 
thickness object in accordance with a preferred aspect of the invention; 
25 Fig. 15 is a flowchart illustrating a method of generating 

phantom landmarks for the phantom thickness map object of Fig. 13; and, 



Fig. 16 is a flowchart illustrating a method of determining breast 
landmarks of the breast thickness object in accordance with an aspect of the 
invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE 
5 INVENTION 

Referring to Fig. 1, there is illustrated in a perspective view, a 
mammography machine 12. The mammography machine 12 includes a 
breast support plate 14, a breast compression plate 18, and an x-ray tube 16. 
In operation, the x-ray tube 16 projects x-rays through the breast compression 

10 plate 18, which is transparent to light and x-rays, through the breast, and 
through the breast support plate 14. The breast compression plate 18 may be 
vertically adjusted to accommodate breasts of different dimensions. The 
breast support plate 14 includes a detector (shown in Fig. 4) that is sensitive 
to the x-rays. Variation in the density and thickness of the breast will affect the 

15 x-rays traveling through the breast. This in turn will affect the image left on the 
detector in the breast support plate 14. These signal variations may then be 
examined for possible malignancies or other conditions. However, to 
determine density, and thus to properly interpret the image, the thickness of 
the breast must be known at all points. 

20 Referring to Fig. 4, a breast that is compressed between the 

breast support plate 14 and the breast compression plate 18 is shown in a 
schematic view. The breast 13 is of a thickness T in centimeters. X-rays 
originating from an x-ray tube 16 project through the breast compression plate 
18, any empty space surrounding the breast, the breast, and breast support 

25 plate 14 to impinge on a detector 20 underneath the breast support plate. 
When they impinge on the detector 20, the x-rays 15 contain information 
about the thickness and composition of that portion of the breast through 
which they have passed. However, the x-rays will also have been affected by 
the empty space between the breast support plate 14 and breast compression 

30 plate 18 through which they have passed. At some points, of course, where 
the breast is in contact with both the breast support plate 14 and breast 



-6- 

compression plate 18, there will be no empty space. However, at other points, 
the curvature of the breast creates a space between the breast compression 
plate 18 and the breast support plate 14 that is not occupied by the breast. If 
the thickness of the breast is known, then this information can be taken into 
5 account when interpreting the x-ray data on the detector 20. 

Referring to Fig. 7, there is illustrated a phantom thickness map 
object 22 generated using three-dimensional modeling software in 
accordance with an aspect of the present invention. This phantom thickness 
map object 22 is generated using a breast phantom 24 constructed of poly- 

10 methyl-methacrylate (PMMA) shown in Fig. 2. This phantom breast 24 is first 
imaged by the mammography machine 12 to obtain a phantom mammogram. 
As the composition of the phantom mammogram image is uniform and known, 
the intensity of the x-rays transmitted through the phantom breast 24 will vary 
based on the variation in the thickness of the phantom breast 24. 

15 Referring to Fig. 3, there is illustrated in a perspective view, a 

three-dimensional triangular phantom 26 in accordance with an aspect of the 
invention. This triangular phantom 26 contains slabs of PMMA 26a, as well as 
plastic layers 26b and 26c simulating 30% and 50% fibroglandular tissue 
respectively (from left to right - PMMA, 30%, 50%, PMMA). This three- 

20 dimensional triangular phantom 26 is then subjected to a mammogram to 
generate a set of image data. Again, this image data will vary only with the 
thickness of the triangular phantom 26. However, the thickness of the 
triangular phantom 26 will be known at any point. Thus, the set of image data 
for the triangular phantom 26 can be used to correlate the thicknesses of the 

25 triangular phantom 26, with particular points in the phantom mammogram 
having the same intensity of x-ray transmitted, and therefore being of the 
same thickness. 

From the x-ray profile along the PMMA triangular phantom 
image from top seven centimeters to base less than one millimeter, the 
30 position along the wedge (i.e. the thickness) is determined from the logarithm 
of the image pixel signal by interpolation using a second-degree polynomial 



fit. This fit is plotted as line 501 on the graph of Fig. 5. The polynomial function 
represented by line 501 of Fig. 5 allows direct conversion from logarithmic 
gray pixel value to thickness value. 

In a similar way, second-degree polynomial functions are found 
5 for 30% fibroglandular tissue and for 50% fibroglandular tissue. The second- 
degree polynomial function for 30% fibroglandular tissue is plotted as line 503 
in Fig. 5, and the second-degree polynomial function for 50% fibroglandular 
tissue is plotted as line 502 in Fig. 5. A second-degree polynomial function for 
100% fibroglandular tissue was obtained by mirroring the 30% polynomial 
10 function around the 50% polynomial function, and is represented as line 500 
in Fig. 5. Any percentage glandular composition can be verified by using 
slabs of known thickness and composition. 

The phantom thickness map 22 and polynomial functions 500, 
501, 502 and 503 can then be used to compute the thickness and density 

15 map of a particular digital mammogram. First, this will require the phantom 
thickness map 22 to be rescaled to the size of the digital mammogram, and 
will require the thickness values of the phantom thickness map 22 to be 
normalized to the thickness readout of the mammographic system. Next, the 
phantom thickness map is overlaid on a digital mammogram image using a 

20 point-based elastic warping method, which is efficient at recovering local 
deformations (see F. Bookstein, Thin-Plate Splines and the Decomposition of 
Deformations, IEEE Transactions Pattern Analysis and Machine Intelligence, 
11, pp. 567-585, 1989). With this technique, special care is needed in the 
selection of landmarks. Two different sets of landmarks are chosen, both in 

25 the phantom thickness map 22 and in the digital mammogram. 

The phantom thickness object of Fig. 7 is defined and 
determined as the object with thickness values larger than zero. Referring to 
Fig. 6, there is illustrated an intensity histogram 29 of a digital mammogram. 
This intensity histogram is bimodal. The breast thickness object 30 of Fig. 8 is 
30 automatically generated from the histogram using a threshold value 32 shown 
in Fig. 6. This threshold value 32 is at the middle point of the valley between 
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the two modes in the histogram. The boundaries of both the phantom 
thickness object of Fig. 7 and the breast thickness object of Fig. 8 are found 
by employing a morphological removing operation. In the binary images of 
Fig. 7 and Fig. 8, a pixel is set to zero (black) if all of its four-connected 
5 neighbours are one (white), thus leaving only the boundary pixels on. 

Having generated the phantom thickness object 22 of Fig. 7, 
and the breast thickness object 30 of Fig. 8, it is possible to select landmarks. 
Referring to Figs. 8 and 9 respectively, the phantom thickness object 22 and 
breast thickness object 30 are shown divided into segments. In the case of 

10 the phantom thickness object 22 of Fig. 9, these segments are defined by a 
series of radial lines 34 extending from the center 32 of the phantom 
thickness object 22 to the outer edge of the phantom thickness object 22. 
Each of these radial lines intersects a first phantom boundary line 38 marking 
the outer edge of the phantom thickness object 22. Together, these 

15 intersection points provide a first set of phantom landmarks 36. Similarly, in 
the case of the breast thickness object 30 of Fig. 10, the segments are 
defined by a series of radial lines 44 extending from the center 42 of the 
breast thickness object 30 to the outer edge of the breast thickness object 30. 
Each of these radial lines intersects a first breast boundary line 48 marking 

20 the outer edge of the breast thickness object 30. Together, these intersection 
points provide a first set of breast landmarks 50. 

Referring to Fig. 9, a second phantom boundary line inside the 
first phantom boundary line is shown. This boundary line represents the point 
at which the breast phantom 24 is no longer in contact with the breast 

25 compression plate 18. This point is selected from a phantom thickness profile 
60 of Fig. 11. Each of the radial lines 34 of Fig. 9 has an associated thickness 
profile such as the thickness profile 60 of Fig. 11. A line 62 is drawn 
connecting the first point 63 and last point 64 of the thickness profile 60. A 
point 66 on the thickness profile 60 is then selected to be a maximum 

30 distance from the line 62. This point 66 is substantially at the point where the 
radial line 34 of the phantom ceases being in contact with the breast 



-9- 

compression plate 18. Together, the points 66 selected for all of the radial 
lines 34, generate the second boundary line 40. 

A breast thickness profile 70 is plotted for each of the radial lines 
44 of the segmented breast thickness object 30 of Fig. 10. The logarithmic 
5 profile values are converted to thickness using the polynomial function for 
50% dense material. A thickness profile 70 of one such radial line 44 is shown 
in Fig. 12. Unlike the thickness profile 60 of the breast phantom 24, the 
thickness of an actual breast is not uniform over a first interval, but instead 
increases before decreasing. Similar to the case of Fig. 11, a line 72 

10 connecting the first point 73 on the profile 70 with the last point 74 on the 
profile 70 is drawn. Then, a point 76 is selected to be a maximum distance 
from the line 72. These points 76 for all of the radial lines 44 are then plotted 
as points 52 on the segmented breast thickness object 30 of Fig. 10, and are 
connected to provide the second breast boundary line 46. Unlike the second 

15 phantom boundary line 40 of the phantom object of Fig. 9, the second 
boundary line 46 of the breast thickness object of Fig. 10 is irregular, 
reflecting variation in the composition and compressibility of the breast. 

The minimum thickness values for thickness on the outer edge 
of the breast thickness object are computed using the polynomial function for 

20 100% dense material, to convert logarithmical grey pixel values to thickness. 
The polynomial function for 100% dense material is selected due to the layer 
of skin surrounding the breast. A corrected warped thickness map is then 
computed by cropping the radial lines 44 and cropping the map generally, at 
the minimum thickness value given by the 100% conversion function. Next, 

25 the cropped profile is approximated by a linear combination of two 
exponentials using a non-linear least squares logarithm. 

Referring to Figs. 13 and 14, there is illustrated in flowcharts a 
method of generating a breast thickness map in accordance with an aspect of 
the present invention. In step 80 of the flowchart of Fig. 13, a phantom 
30 mammogram is obtained by imaging a breast phantom 24. This phantom 
mammogram contains a series of profiles of the breast phantom along 
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different planes, reflecting the difference in thickness of the breast phantom at 
these different planes. The image data for different thicknesses is then 
generated by imaging a three-dimensional triangular phantom 26. This 
triangular phantom contains slabs of PMMA, as well as plastics simulating 30 
5 and 50% of fibroglandular tissue. By imaging this three-dimensional triangular 
phantom 26, image data for known and different thicknesses are generated. 
This information can then be combined with the information provided by step 
80, to determine a phantom thickness map object 22 in step 84. Then, in 
steps 86 and 88, a first set of phantom landmarks, and a second set of 

10 phantom landmarks are determined. 

Referring to Fig. 14, a digital mammogram of an actual breast is 
obtained in step 90. Then, in steps 92 and 94 respectively, a first set of breast 
landmarks, and a second set of breast landmarks are defined. In step 96, the 
phantom thickness map is rescaled to the size of the digital mammogram, and 

15 in step 98, the phantom thickness map object is normalized by normalizing its 
thickness size to the thickness readout of the mammography system. In step 
100, the phantom thickness map object is overlaid on the digital mammogram 
using a point-based elastic warping method, which is efficient at recovering 
local deformations. 

20 Referring to Fig. 15, there is illustrated in a flowchart a method 

of selecting a first set of phantom landmarks and a second set of phantom 
landmarks in accordance with an aspect of the present invention. In step 110, 
a series of radial lines extending from the center of the phantom thickness 
object 22 to its outer edge are generated. In step 112, a first set of phantom 

25 landmarks are determined by taking the intersection of these radial lines with 
the outer edge of the phantom thickness object. Then, in step 114, a 
secondary boundary of the phantom thickness object 22 is determined. This 
secondary boundary of the phantom thickness object is defined by the points 
at which the phantom thickness object moves from being in contact with the 

30 breast compression plate, to not being in contact with the breast compression 
plate. Then, in step 116, a second set of phantom landmarks is determined. 
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This second set of phantom landmarks is determined by taking the 
intersection of the radial lines generated in step 110 with the secondary 
boundary generated in step 114. 

Referring to Fig. 16, there is illustrated in a flowchart a method 

5 of selecting a first set of breast landmarks and a second set of breast 
landmarks in accordance with an aspect of the present invention. In step 120, 
a series of radial lines extending from the center of the breast image to its 
outer edge are generated. In step 122, a first set of breast landmarks are 
determined by taking the intersection of these radial lines with the outer edge 

10 of the breast image. Then, in step 124, a secondary boundary of the breast 
image is determined. This secondary boundary of the breast image is defined 
by the points at which the breast changes from being in contact with the 
breast compression plate, to not being in contact with the breast compression 
plate. Then, in step 126, a second set of breast landmarks is determined. This 

15 second set of breast landmarks is determined by taking the intersection of the 
radial lines generated in step 120 with the secondary boundary generated in 
step 124. 

According to a preferred aspect of the present invention, step 
100 of the flowchart of Fig. 14 is executed by applying the point-based elastic 
20 warping method to warp the first set of phantom landmarks into the first set of 
breast landmarks, and to warp the second set of phantom landmarks into the 
second set of breast landmarks. 

Other variations and modifications of the invention are possible. 
For example, phantom thickness objects may be generated in other ways by, 
25 say, for example, assembling an average breast from a series of 
mammograms for different women, or by selecting a stored breast thickness 
object that most closely matches the shape and dimensions of the breast 
being imaged from a library of previously obtained breast thickness objects. 
Further, other techniques may be applied to overlay the phantom thickness 
30 map on the breast thickness object. All such modifications or variations are 



believed to be within the sphere and scope of the invention as defined by the 
claims appended hereto. 



